UMHEP-427 



Systematics of soft final state interactions in B decay 

John F. Donoghue, Eugene Golowich, Alexey A. Petrov and Joao M. Scares 
Department of Physics and Astronomy, University of Massachusetts, Amherst MA 01003 USA 

Abstract 

By using very general and well established features of soft strong interac- 
tions we show, contrary to conventional expectations, that (i) soft final state 
interactions (FSI) do not disappear for large m^, (h) inelastic rescattering is 
expected to be the main source of soft FSI phases, and (iii) flavor off-diagonal 
FSI are suppressed by a power of m^, but are quite likely to be significant at 
niB — 5 GeV. We briefly discuss the infiuence of these interactions on tests 
of CP-violation and on theoretical calculations of weak decays. 
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It is notoriously difficult to say anything useful about final state interactions in weak 
decays. Although the final state interactions are not themselves of fundamental interest, 
they are important for some truly interesting aspects of B decay. For example, many signals 
of direct CP violation in B transitions require final state phases as well as CP-violating 
phases if the CP-odd asymmetry is to be nonzero. In this paper we shall derive some 
general properties of soft final state interactions and describe the implications for theory 
and phenomenology. 

The scattering of hadrons at high energies exhibits a two-component structure of 'soft' 
and 'hard' scattering. Soft scattering is that which occurs primarily in the forward direction. 
The transverse momentum is limited, having a distribution which falls exponentially with 
a scale of order 0.5 GeV. At higher transverse momentum ultimately, one encounters the 
region of hard scattering, which falls only as a power of the transverse momentum. Collisions 
involving hard scattering are interpreted as interactions between pointlike constituents of the 
hadrons, the quarks and gluons of QCD. These are calculable in QCD perturbation theory 
and are found to be in good quantitative agreement with experiment. Hard scattering is, 
however, only a very small portion of the total hadronic cross section. The much larger 
soft component at low values of transverse momentum is by far the dominant contribution 
to high energy scattering. Although soft hadronic interactions are generally not calculable 
from ffist principles, there is available a wealth of experimental studies ^ and accurate high 
energy phenomenology ^ on which to base our study. 

The modern approach to B physics employs as an organizing principle the fact that the 
B mass is very large compared to the QCD scale. In the context of soft FSI in B decays, 
it suggests the question — what is the leading order behavior of soft final state phases in 
the ttlb —>■ oo limit? The common perception among researchers is that they should become 
less and less important as the mass of the decaying quark becomes heavier. This is because, 
roughly speaking, 'the final state particles emerge at such high momenta that they do not 
have a chance to rescatter'. Such an expectation is, however, false because soft scattering 
actually grows with energy. As an example of this important energy dependence, we shall 
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vtscTs^f^^^ (/|r|5)-(/|rt|B)] =l^(/|rt|j)(j|r|B) . (i) 



demonstrate below that the imaginary part of the forward elastic amplitude has an s^^'^ 
{rj ~ 0.08) dependence, and as a consequence, the elastic final state interaction is roughly 
constant as a function of rriB- We shall then use this observation as the starting point for 
a more general exploration of the systematics of FSI for large rriB- The inevitability of our 
conclusions will be seen to follow rather directly from well established aspects of strong 
interaction phenomenology. 

Final state interactions in B decay involve the rescattering of physical final state particles. 
Unitarity of the 5-matrix, iS^iS = 1, implies that the T- matrix, iS = 1 + iT, obeys 
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Of interest are all physical intermediate states which can scatter into the final state /. Among 
all these, however, we shall first concentrate on just the elastic channel and demonstrate 
that elastic rescattering does not disappear in the limit of large m^.Q The elastic channel is 
especially convenient for our discussion because we can use the optical theorem to rigorously 
connect it to known physics. The optical theorem relates the forward invariant amplitude 
to the total cross section, 

Im Mf^f{s, t = 0) = 2k^/saf^a.n ^ saf^^ii , (2) 

where s is the squared center-of-mass energy and t is the squared momentum transfer. 

The asymptotic total cross sections are known experimentally to rise slowly with energy. 
All known cross sections can be parameterized by fits of the form @ 

g \ 0.08 / g \ -0.56 



where Sq = 0{1) GeV is a typical hadronic scale. Thus, the imaginary part of the forward 
elastic scattering amplitude rises asymptotically as This growth with s is counterin- 



^We stress that we are not suggesting the elastic channel to be the dominant contribution to soft 
rescattering. Our analysis leads to quite the opposite conclusion, that it is the inelastic channels 
which are most important. 



tuitive in that it cannot be generated by a perturbative mechanism at any finite order. In 
particular, calculations based on the quark model or perturbative QCD would completely 
miss this feature. 

In order to arrive most simply at our goal, let us first consider only this imaginary part, 
and build in the known exponential fall-off of the elastic cross section in t (recalling that t 
is negative) by writing 



s 



1.08 
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ilm Mf^fis, t) ~ z/3o — e''* . (4) 

KsqJ 

It is then an easy task to calculate the contribution of the imaginary part of the elastic 
amplitude to the unitarity relation for a final state / = a + b with kinematics p'^+Pb = Pa+Pb 
and s = {pa + PbY, and we find 

where we have used t = {pa — p'a)"^ — —s{l — cos^)/2 and have taken s = m\. The 
integration over the angle involving the direction of the intermediate state is seen to introduce 
a suppression factor to the final state interaction of = m']^. This is because the soft final 
state rescattering can take place only if the intermediate state has a transverse momentum 

< 1 GeV with respect to the final particle direction. This would naively suggest a result 
consistent with conventional expectations, i.e. an FSI which falls as m^^. However, the fact 
that the forward scattering amplitude grows with a power of s overcomes this suppression 
and leads to elastic rescattering which does not disappear at large mB- 

In fact, we can make a more detailed estimate of elastic rescattering because the phe- 
nomenology of high energy scattering is well accounted for by Regge theory. % Scattering 
amplitudes are described by the exchanges of Regge trajectories (families of particles of 
differing spin) which lead to elastic amplitudes of the form 

/ Q \ a(t) 

Mf^f=mt)[jJ (6) 
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with ^ = 1 for charge conjugation C = +1 and ^ = i for C = — 1. Each such trajectory is 
described by a straight hne, 



The leading trajectory for high energy scattering is the Pomeron, having C = +1, ao ^ 1.08 
and a' ^ 0.25 GeV~^. Note that since 



the exponential fall-off in t is connected with the slope a' and the effective slope parameter 
b in Eq. thus increases logarithmically with s. Since ao is near unity, the phase of the 
Pomeron-exchange amplitude is seen from Eq. (|^) to be almost purely imaginary. This 
feature has been verified experimentally by interference measurements. There are several 
next-to-leading trajectories, both those with C = —1 (p(770) & c<j(782) trajectories) and 
those with C = +1 (a2(1320) & /2(1270) trajectories). Roughly, these have ao ^ 0.44, 
a' ~ 0.94 GeV~^ and lead collectively to the s~^'^^ dependence in the asymptotic cross 
section of Eq. (^). The prefactor /5(t) in Eq. (|^) also has known regularities. For the 
Pomeron, f3 is very nearly proportional to the number of quarks at each vertex, and carries 
a power law behavior similar to the electromagnetic form factor. Therefore, (3.,^^ in pion-pion 
scattering can be expressed in terms of the analogous proton-proton quantity jSpp as 



The combination of exponential and power law t dependence in a generic Regge amplitude 
gives a unitarity integral no longer having an elementary form. However, the integration 
can still be carried out in terms of Euler functions. Taking s = m\ ~ 25 GeV^, we obtain 
for the Pomeron contribution 



a{t) = ao + Ci't 



(7) 




(8) 




(9) 



Vise Mb 



■TTTT Pomeron 



■TTTT 



(10) 



where we find from our computation. 
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e~0.21 . (11) 

iFiom this numerical result and from the nature of its derivation, we may anticipate that 
additional individual soft FSI will not be vanishingly small. Moreover, other final states 
should have elastic rescattering effects of comparable size. However, of chief significance is 
the weak dependence of e on that we have found — the (m^)° °^ factor in the numerator 
is attenuated by the ln(m^/so) dependence in the effective value of b (compare Eqs. (^,(^). 

The above study of the elastic channel, although instructive, is far from the whole story. 
In fact, it suggests the even more significant result that at high energies FSI phases are 
generated chiefly by inelastic effects. At a physical level, this conclusion is forced on us 
by the fact that the high energy cross section is mostly inelastic. It is also plausible at 
the analytic level, given that the Pomeron elastic amplitude is almost purely imaginary. 
The point is simply this. Our study of elastic rescattering has yielded a T-matrix element 
Tab^ab = 2ie, which directly gives Sab^ab = 1 — 2e. However, the constraint of the iS-matrix 
be unitary can be shown to imply that the off-diagonal elements must be 0{y/e). Since e is 
approximately (9(m^) in powers of and numerically e < 1, the inelastic amplitude must 
also be (9(m^) and of magnitude ^/e > e. There is an alternate argument, utilizing the 
form of the final state unitarity relations, which also shows that inelastic effects are required 
to be present. In the limit of T-invariance for the weak interactions, the discontinuity 
T>isc M.B~^f is a real number (up to irrelevant rephasing invariance of the i?-state). The 
factor of i obtained in the elastic rescattering in Eq. (|10D must be compensated for by the 
inelastic rescattering (this effect is made explicit in the example to follow) in order to make 
the total real. Therefore, the presence of inelastic effects is seen to be necessary. 

Analysis of the final-state unitarity relations in their most general form, 

Vtsc Mb^h = \Y. Ms^kTLf, , (12) 

k 

is quite complicated due to the many contributing intermediate states present at the B 
mass. However, it is possible to illustrate the systematics of inelastic scattering by means 
of a simple two-channel model. This pedagogic example involves a two-body final state /i 
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undergoing elastic scattering and a final state /2 which is meant to represent 'everything 
else'. We assume that the elastic amplitude is purely imaginary. Thus, the scattering can 

be described in the one-parameter form 

( cos 29 2 sin 2^ \ (2ism^e sin 20 \ 

S=[ , T=[ , (13) 

\ism2e cos 29 J \ sin 20 2ism^eJ 

where, from our elastic-rescattering calculation, we identify sin^ 9 = e. The unitarity rela- 
tions become 

1 
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Disc M B^fi = -i sin^ 9M + 2 ^^"^ ^^^2 ' 



Vise M B~*f2 = ^ sin 26'A< B^h ~ ^ sin^ 9M B^f2 (14) 

If, in the limit 0^0, the decay amplitudes become the real numbers }A\ and these 
equations are solved by 

Mb-^h = cos 0A^? + i sin 9Ml , -Ms^/a = cos 9Ml + i sin 0A^° . (15) 

As a check, we can insert these solutions back into Eq. (p^. Upon doing so and bracketing 
contributions from Mb^/i and M.B^f2 separately, we find 
1 



Vise Mb^/i = - 



(16) 



The first of the four terms comes from the elastic channel /i and is seen to be cancelled by 
the final term, which arises from the inelastic channel /2. The third term is dominant, being 
0{^/e), and comes from the inelastic channel. 

In this example, we have seen that the phase is given by the inelastic scattering with a 
result of order 

ItuMb^^ (17) 
UeMB^f ^ Ml ■ ^ ' 

Clearly, for physical B decay, we no longer have a simple one-parameter S matrix. However, 
the main feature of the above result is expected to remain — that inelastic channels cannot 
vanish because they are required to make the discontinuity real and that the phase is sys- 
tematically of order ^/e from these channels. Of course, with many channels, cancellations 



or enhancements are possible for the sum of many contributions. However the generic expec- 
tation remains — that inelastic soft final-state-rescattering arising from Pomeron exchange 
will generate a phase which does not vanish in the large ttib limit. 

What about nonleading effects? It is not hard to see that these may be significant at 
the physical values of m^. For example, the fit to the pp total cross section is 

. X r / S \ 0.08 / S \ -0.56i , , , ^ 

a{pp) = 22.7( — ) + 140(— ) {mb) (18) 

•So ^0 

with So = 1 GeV^. At s = (5.2 GeV)^, the nonleading coefficient is a factor of six larger 
that the leading effect, effectively compensating for the s~^-^^ = m'^'^^ suppression. The 
subleading terms are then comparable in the elastic forward pp scattering amplitude. The 
slope of the p trajectory and hence the experimental fall-off with t, is larger than that of 
the Pomeron by a factor of nearly four, and thus this moderates the integrated rescattering 
effects. If we estimate the (5 coefficient of the p trajectory in tttt by relating it to pp via a 
factor of ~ 4/3^^ and then perform the integration over the intermediate state momentum 
we find 



Disc Mb- 



= iepMB^nn , (19) 

p-traj 



with Cp ~ 0.11 — 0.05 i. It is likely that the /2(1270) trajectory could be somewhat larger, 
as it is in pp and irp scattering. 

Final state phases can contribute to weak decay phenomenology in a variety of ways. 
Here, we briefly consider two of these, isospin sum rules and CP- violating asymmetries. A 
simple example of an isospin sum rule is the following relation between B ^ tttt decay 
amplitudes, 

2\/2 

A^+_ - Moo = -^M-o , (20) 

where M.+- = M.{B^ — > tt+tt"), etc. Measurement of the magnitude of each amplitude via 
the partial decay rate allows one to test the sum rule. Noting that the tttt final state in B 
decay occurs in the isospin states / = 0, 2, one can solve for the difference in phase angles. 



cos((5o - ^2) = 7 , (21) 

4 |A^_o|a/9|A^+-|2 + 9\Moo\^ - 4|^_oP 



At a theoretical level, one sees that the leading Pomeron effect does not contribute to these 
isospin sum rules since Pomeron exchange is identical for each vrV-' final state and thus 
generates only a common overall phase. Thus, the phases measured in isospin sum rules are 
technically subleading, of order m^^'^^. 

CP-violating asymmetries involve comparisions of B ^ f and B ^ f . In order to 
be nonzero, these require two different pathways to reach the final state /, and these two 
paths must involve different CP-violating weak phases and different strong phases. The 
leading Pomeron phases can contribute to such asymmetries if the other conditions are 
met. Because the strong phase is generated by inelastic channels, the relevant pathways 
would involve B f directly or B —>■ 'multibody' followed by the inelastic rescattering, 
'multibody' —* f. Depending on the dynamics of weak decay matrix elements, these may 
pick up different weak phases. As an example, consider the final state / = i^~7r°, which can 
be generated either by a standard W exchange or by the penguin diagram, involving different 
weak phases, For the strong rescattering, we must also consider a channel to which K~7r^ 
scatter inelastically, which we call Knir (although one can generate this asymmetry by a 
hard rescattering DgD K^tt'^, we are concentrating here on the soft physics). The W- 
exchange and penguin amplitudes will contribute with different weight to Kir and Knir, so 
that in the absence of final state interactions we expect 

M{B~ K~7i) = \ Ai\e''^' = AXe''^^" + A\e'^^ 
M{B- K-rnr) = |A„|e^'^" = A'^e'^^ + A^e'^^ (22) 

with 01 7^ (f)n- If we now model the strong rescattering by the two channel model described 
above, we have for B and B decays 

M{B- K'tt) = \Ai\e''^'+i^ \Arr\e''''" 
M{B+ K+uTc) = \Ai\e-'^' + |A„|e-*'^" (23) 
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This leads to a CP-violating decay rate asymmetry 

r(B- ^ X-7r°) - r(B+ ^ X+7r°) ~ v^|Ai||A„| sin((/)„ - 0i) (24) 

While this effect will be very difficult to calculate, we see that inelastic final state interactions 
can contribute to CP- violating asymmetries at leading order in m^. 

The results obtained in this paper must also be accounted for in any theoretical cal- 
culation of weak decay amplitudes. For large niB, there is the hope that one can directly 
calculate the weak matrix elements through variants of the factorization hypothesis or by 
perturbative QCD. Final state interactions will impose limits on the accuracy of such meth- 
ods, as no existing technique includes the effect of inelastic scattering. There must exist, 
in every valid theoretical calculation, a region of the parameter space where the nonpertur- 
bative Regge physics is manifest. Arguments based on local quark-hadron duality do not 
account for these effects of soft physics because the growth of the scattering amphtude with 
s (for both the leading and first nonleading trajectories) cannot be seen in perturbative cal- 
culations. It remains an intriguing possibility that the assumption of quark-hadron duality 
can be questioned in other aspects of S-decay as well. At any rate, for final state interac- 
tion studies, one may only hope that the perturbative/calculable physics is larger then the 
difficult nonperturbative contributions discussed in this paper. 

To conclude, we have argued that the general features of soft scattering have forced 
upon us some suprising conclusions regarding final state interactions. Most importantly, the 
growth of forward scattering with s, as required by the optical theorem and cross section 
data, indicates that soft scattering does not decrease for large m^. The structure of the 
elastic rescattering via the Pomeron also requires that inelastic processes are the leading 
sources of strong phases. These systematics can be important for the phenomenology of B 
decays. 
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